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We investigate the non-adiabatic dynamics of spinning black hole binaries by using an analytical 
Hamiltonian completed with a radiation-reaction force, containing spin couplings, which matches 
the known rates of energy and angular momentum losses on quasi-circular orbits. We consider both 
a straightforward post-Newtonian-expanded Hamiltonian (including spin-dependent terms), and a 
version of the resummed post-Newtonian Hamiltonian defined by the Effective One-Body approach. 
We focus on the influence of spin terms onto the dynamics and waveforms. We evaluate the energy 
and angular momentum released during the final stage of inspiral and plunge. For an equal-mass 
binary the energy released between 40 Hz and the frequency beyond which our analytical treatment 
becomes unreliable is found to be, when using the more reliable Effective One-Body dynamics: 
0.6%M for anti-aligned maximally spinning black holes, 5%M for aligned maximally spinning black 
hole, and 1.8%M for non-spinning configurations. In confirmation of previous results, we find that, 
for all binaries considered, the dimensionless rotation parameter J/E^ is always smaller than unity 
at the end of the inspiral, so that a Kerr black hole can form right after the inspiral phase. By 
matching a quasi-normal mode ringdown to the last reliable stages of the plunge, we construct 
complete waveforms approximately describing the gravitational wave signal emitted by the entire 
process of coalescence of precessing binaries of spinning black holes. 

PACS numbers: 



I. INTRODUCTION 

An international network of kilometer-scale laser-interferometric gravitational-wave detectors, consisting of the 
Laser-Interferometer Gravitational- wave Observatory (LIGO) Q, of VIRGO |3], of GEO 600 Q and of TAMA 300 
0, has by now begun the science operations. TAMA 300 reached its full design sensitivity in 2001, VIRGO is in its 
commissioning phase and plans to start the first scientific runs by the end of 2005, while LIGO has already completed 
three science runs (two of them in coincidence with GEO 600) with increasing sensitivity and duty cycle. LIGO and 
GEO 600 are expected to reach their full design sensitivity in 2005. 

Binary black holes are among the most promising sources for these detectors. Among black hole binaries, it was 
emphasized in Q that there is a bias towards first detecting mostly aligned spinning binaries with high masses, 
whose last stable spherical orbits are drawn, by spin effects, to larger binding energies, yet due to their high masses 
these energies are still emitted through gravitational waves in the sensitive band of the detectors. Studying in detail 
the waveforms emitted during the last stages of dynamical evolution of such heavy spinning black hole binaries, 
with explicit consideration of the crucial transition between adiabatic inspiral and plunge, is a demanding theoretical 
challenge. The aim of the present paper is to provide a first attack on this problem by using some of the best analytical 
tools currently available to describe the transition from adiabatic inspiral to plunge, and notably the Effective One 
Body (EOB) approach HQ. 

So far, most theoretical and data-analysis studies on precessing binaries of spinnin g bl ack holes assumed adiabatic 
evolution. Thus they were restricted to considering only the inspiral phase H, 0, ITonTTl IT^ IT^ lisl IT^ IT^ [l8l | . 
Actually, even for non-spinning binary configurations, most theoretical studies confined themselves to considering the 
adiabatic inspiral phase, though a lot of effort was spent to improve the accuracy of the phasing during the last stages 
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of the inspiral, see e.g. [T^l20| . 



For heavy black hole binaries, most of the signal-to- noise ratio in the ground detectors will come from the very last 
stages of the inspiral, and from the non-adiabatic transition between inspiral and plunge. It is therefore essential to 
be able to describe, with acceptable accuracy, this non-adiabatic evolution. In Refs. '^,t| a new way of describing the 
dynamics of binary systems was introduced: the Effective One-Body (EOB) approach. The EOB approach uses both 
a specific resummation of the post-Newtonian Hamiltonian, and a resummed version of radiation reaction. This was 
shown to lead to a rather robust formalism, which is likely to provide a reliable description of non-adiabatic effects, 
of the transition between inspiral and plunge, and of the beginning of the plunge. It was also used in Jj to model the 
full merger phase of non-spinning binaries, by matching the natural end of the EOB plunge with the ring-down phase. 
The EOB approach was used in Refs. l2lj| to derive non-adiabatic template waveforms emitted by non- spinning 
black hole binaries. It was shown in that these new EOB templates led to significantly enhanced signal to noise 
ratios in current detectors (mainly because of the inclusion of the plunge signal) . The EOB Hamiltonian was extended 
to the case of spinning black holes in (5| . The analytical predictions made by the EOB method (including spin) were 
found to agree remarkably well |23| with the numerical results obtained by means of the helical Killing vector approach 
[j^l for circular orbits of corotating black holes. Several other studies showed that the EOB method provides phasing 
models which are more reliable and robust than other (adiabatic or non-adiabatic) models 25. .26] . 

The main purpose of this paper is to extend the use of the EOB approach to the case of precessing binaries 
of spinning black holes, both by including spin-dependent terms in radiation-reaction effects, and by studying the 
waveforms generated beyond the adiabatic approximation, i.e. taking into account the transition between inspiral 
and plunge, and the plunge itself. Let us emphasize again that the EOB approach has the advantage of providing 
an analytical description of the transition from inspiral to plunge. Recently, some attempts have been made to 
tackle, by means of 3D numerical simulations (combined with a p erturbative approach), the gravitational radiation 
emitted by a very tight black hole binary both in non-spinning '275 j ^^'^ moderately spinning, but non-precessing 



configurations. These three-dimensional (3D) simulations concluded to the emission of a significantly larger amount 
of energy in the form of gravitational radiation than what we shall find from our analytical, EOB approach. It should 
be mentioned in this respect that the energy released in the form of gravitational waves depends very much on the 
choice of initial data, and that the amount by which the initial data chosen in |27t .28.] differ from the physically 
correct "no-incoming-radiation" data is unclear. This crucial issue will be further discussed below. 

Recent simulations 18] based on population synthesis codes predict that ~ 50% — 80% of neutron star-black hole 
(NS-BH) binaries in the Galactic field may have tilt angles (i.e., the angle between the black hole (BH) spin and the 
orbital angular momentum) between and 40° and 10 — 20% of NS-BH binaries between 40° and 50°. By studying 
the formation of close compact binaries and the misalignment angle that can occur after the second core-collapse 
event, Kalogera (2^ predicted that the majority of BH-BH binaries in Galactic binaries may have a tilt angle smaller 
than 30°. All these results assume that the misalignement is entirely due to the recoil velocity ("kick") imparted to 
the NS (or the smaller BH in the binary) at birth by the core-collapse. However, the spin properties of NS-BH and 
BH-BH binaries in dense environment and centers of globular clusters could be very different than in the Galactic field. 
Considering the low event rates, ~ 1 per 2 years, of binary coalescences in first generation of ground-based detectors, 
it is worthwhile to adopt a conservative point of view and investigate waveforms for generic spin configurations. Little 
is known about the magnitude of the spin of NSs and BHs. From the observed pulsars the dimensionless rotation 
parameter ons takes values in the range 0.005 — 0.02. Our analysis will focus on BH-BH binaries and we shall consider 
arbitrary spins: < obh < 1- 

For completeness, we investigate the two-body dynamics by adopting two approaches: the straightforward post- 
Newtonian (PN)-expanded Hamiltonian l30Ll3ll] and the PN-resummed Hamiltonian a la EOB 0, j2j, ^ [henceforth 
referred to simply as the EOB-Hamiltonian] . For simplicity, instead of using the Kerr-deformed EOB-Hamiltonian 
derived by Damour in Ref. 5], we shall use as Hamiltonian in this paper the sum of the purely orbital (Schwarzschild- 
deformed) EOB-Hamiltonian i21j and of the separate contributions due to spin-orbit and spin-spin effects. [Note 
that, among the spin-spin terms, there are the terms due to monopole-quadrupole interactions js^], Q.] 

The paper is organized as follows. Section ITTl summarizes the main results of the conservative part of the two-body 
dynamics in the Hamiltonian formalism, and contains formula for the PN-expanded and EOB Hamiltonians up to 
3PN order. In Sec. IIIII we augment the dynamics with radiation-reaction effects. We derive the radiation-reaction 
force which includes spin effects and matches known rates of energy and angular momentum losses on quasi-circular 
orbits. [Our result agrees with the recent results of Will |3^] that appeared after we had completed our work.] In 
Sec. IIVI we define the two-body approximants and discuss the initial conditions used to evolve the precessing two- 
body dynamics; we compare the (lowest-order) waveforms obtained using PN-expanded and (PN-resummed) EOB 
dynamics by computing the overlaps between these two types of waveforms. In view of the greater a priori reliability 
of the EOB approach, we use them (and only them) to estimate the energy and angular momentum released during 
the last stages of evolution. Section IVlII contains our main conclusions. 

We leave to future work a thorough application of our results to data analysis purposes. 
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II. CONSERVATIVE HAMILTONIAN INCLUDING SPIN-ORBIT AND SPIN-SPIN EFFECTS 
A. Orbital third-post-Newtonian expanded Hamiltonian 

The purely orbital (non-spinning)third-post-Newtoman Hamiltonian H'^ = H'^{X, P) (in the center of mass frame, 
and after subtraction of the total rest-mass term) was derived in Ref. 31] (completed by Refs. |3^|35|). In scaled 
variables, and written as a straigthforward PN-expansion, it reads (see Ref. 21]): 



i7°(x,P)Exp = f^H°{q,p) = ^ #Ncwt(q,p) + i^iPN(g,p) + H2PN{q,p) + i?3PN(q,p) 



(2.1) 



where fj, = mi 1712/ M, M = mi + TO2 and (q, p) denote the canonical variables p = Pi//i = — P2//i, and q = X/il/ — 
(Xi — X2)/M, where Xi and X2 are the positions of the black hole centers of mass in quasi-Cartesian Arnowitt-Deser- 
Misner (ADM) coordinates. The Newtonian term and the PN contributions read (denoting rj = /i/M ~ mim2/M^): 
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where the scalars q and p are the (coordinate) lengths of the two vectors q and p; and the vector n is just q/g. 
B. Orbital third-post-Newtonian efFective-one-body Hamiltonian 



(2.5) 



(2.6) 



As was emphasized in previous work (see e.g. |36l l37| '). and as we shall confirm below, the non-resummed PN- 
expanded Hamiltonian (or the non-resummed PN-expanded equations of motion) do not lead to a reliable description 
of the evolution near the last stable circular orbit, nor, a fortiori during the transition between inspiral and plunge. On 
the other hand, it was argued in [1,13,0] that the EOB approach defines a specific resummation of the PN-expanded 
Hamiltonian which leads to a much more reliable description of the dynamical evolution near the last stable circular 
orbit, and of the transition between inspiral and plunge. 

The explicit expression of the purely orbital, EOB-Hamiltonian is jQ]: 



where H^s is given by 0, |^ : 

iIeff(X',P') = M^cff(q',p') 



M. 



(2.7) 
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with q' and p' being the reduced canonical variables obtained by rescaling X' and P' by M and /i, respectively, 
n' = q'/q' where we set q' = |q'|. The coefficients 21,2:2 and 23 are arbitrary, subject to the constraint 

821 +422 + 323 = 6(4 - 3??) 77. (2.9) 

Setting (as in Ref. '25^) 21 = r]Zi, 22 = 7722 and 23 = 7723, the coefficients A{q') and D{q') in Eq. (|2.8|) read: 

'94 41 



q' q'-^ 
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21 



D{q') = 1 - ^ + [721 + 22 + (377 - 26)] ^ 



(2.10) 
(2.11) 



As done in Ref. we restrict ourselves to the case 21 = 22 = and improve the behavior ^ by replacing the 
"effective potential" A{q') with the Fade approximants 
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C. Adding spin couplings 

There are several ways of including spin effects in the Hamiltonian dynamics. When considering the PN-expanded 
form of the orbital Hamiltonian H^^p, it is natural to add the spin-dependent terms as further additional contributions: 
Htot = iJ^xp + iJ^P^^ . On the other hand, when considering the EOB-resummed form of the Hamiltonian ^eoB' 
it has been argued in Ref. Q that it is probably better to incorporate most of the spin effects within a suitably 
generalized (a la Kerr) EOB-Hamiltonian, whose explicit form will be found in |5l|. In the present work, we shall, for 
technical simplicity, adopt a uniform way of including spin effects. Namely, we shall simply include them a linearly 
added contributions to the basic (PN-expanded or EOB-resummed) orbital Hamiltonian. We shall verify below that, 
in the EOB-resummed case, the two different ways (a la (sl, or as in the following equation) of incorporating spin 
effects lead to very similar physical effects. 

Finally, the explicit forms we shall use of the "spinning" Hamiltonian read: 

Hexp{X, P, Si, S2) = H°^p{X, P) + Hso{X, P, Si, S2) + Hss{X, P, Si, S2) , (2.15) 

or 

Heob{X, P, Si, S2) = H^oBiX, P) + Hso{X, P, Si, S2) + Hss{X, P, Si, S2) , (2.16) 



^ As shown in |21| . the use of the straightforward PN-expanded, 3Pt'(- accurate "effective potential" A{q') does not lead to a well-defined 
last stable circular orbit (contrary to what happens in the 2PN-accurate case |^). This is due to the rather large value of the 3PN 
coefficient ^ - ^tt^ ~ 18.688 entering the PN expansion of A{q'). Replacing the PN-expanded form of A(g') by a Pade approximant 
cures this problem. [See also l23l for a comparison of the physical consequences of various possible resummations of A{q').] 
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where |30|, IM |39| : 

Hss - Hs,s, + Hs,s, + Hs,s, = ^ [3(^o ■ N)iSo ■ N) - {So ■ S^)] , (2.18) 
So^(l+"-^)s,+ (l+'-^)s,, (2.19) 

Hs,s, = ^ [3(5i • N){S2 ■ N) - (5i • S2)] , (2.20) 
Hs.s, + Hs^s, = [3(Si • Ar)(Si ■ AT) - {S^ ■ S,)] — + ^ [3(^2 • N)iS2 ■ N) - [82 ■ S^)] ^ (2.21) 

with N — X / R, R= \X\ and L = X x P. The spin-spin term -ffsiS2 "^^-^ derived in Ref. [s^, while the spin-spin 
terms Hg-^Si, Hg^s2 which are vahd only for a BH binary were derived in Ref. [see discussion around Eqs. (2.51)- 
(2.55) in Ref. 5] and also Ref. ^3]- They originate from the interaction of the monopole m2 with the spin-induced 
quadrupole moment of the spinning black hole of mass mi and viceversa. The spin-induced quadrupole moment of a 
NS depends on the equation of state. So, if we were applying our approch to NS binaries we could take into account 
the monopole-quadrupole interaction by multiplying Hs-^Sn HS2S2 by some equation-of-state-dependent coefhcient 7 
[see Ref. |,33j]. 



D. Equations of motion and conserved quantities 

The time evolution of any dynamical quantity f{X, P . Si, S2) is given by |^: 



j^f{X,P,Si,S2) = {f,H}, (2.22) 



where with {..., ...} we indicated the Poisson brackets {X^, Pj} = S'j. The Hamilton equations of motion are: 

dX _ dH dP _ dH 

HT^^'dP' ~dt~^'dX' 

The equations of motion for the spins are easily derived as 0, E3| '■ 



with 



(2.23) 



J n rr 

-Si = {Si,H} = — xSi=^ixSi, (2.24) 

-S2^{S2,H} = — X S2 = n2X S2, (2.25) 
dt 0S2 



Note that a consequence of the above spin-evolution equations is the constancy of the lengths of the spin vectors: 
Sl = est., ^2 = est. 

When using the EOB Hamiltonian Eq. H2.16|l we should in principle consider the canonical transformation between 
(X,P) and {X',P') which is explicitly given as a PN expansion in Refs. 0, |^. However, since the Hamilton 
equations are valid in any canonical coordinate system, when we evolve the EOB dynamics we write the Hamilton 
equations in terms of {X',P') and for the spinning part we neglect the differences between {X',P') and {X,P) 
which are of higher PN order. When in the following sections we compare the results between PN-expandcd and 
PN-resummed Hamiltonians, we will always compare quantities which are gauge invariant to lowest order. 

One of the advantages of using an Hamiltonian formalism is that one can immediately derive from the fundamental 
symmetries of the relative dynamics (time translation and spatial rotations) two exact conservation laws: that of the 
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total energy E = H, and that of the total angular momentum J = L + Si + 82- Using Eqs. H2.24(l . 12.25|l it is easy 
to check the conservation of the total energy, 

dH _ OH dX dH dP dH dSi dH dS2 

IT ^ 'dX ^ W IT ^ dSl ~dr ^ dS^ ~dr ^ ' -* 

Similarly, one easily checks the conservation of the total angular momentum J = L + Si + 82- Note the remarkable 
fact that the orbital contribution to J is exactly given, at any PN-order, by the "Newtonian-looking" expression 
L = X X P. This is contrary to what happens when working within a Lagrangian formalism, where the expression 
for the conserved orbital angular momentum gets modified at each PN order: L = fiX x t> + 0(c^^). Here, all the 
needed PN contributions arc included in the Hamiltonian H (and thereby imply that V = dH/dP is a complicated 
function of V). 

E. Spin-orbit interaction and "spherical orbits" 

For most of the dynamical evolution, the spin-spin terms are much smaller than the spin-orbit ones. If we restrict 
ourselves to spin-orbit interactions, the equation of motion for the orbital angular momentum reads 

^ = {L,i/} = ASeff XL. (2.29) 

A useful consequence of this approximate evolution law is that is conserved. Under the same approximation, the 
total spin 8 = 81 + 82 satisfies the equation: 

^^-^^e.xX. (2.30) 

The above (approximate) evolution equations exhibit clearly the (exact) conservation of the total angular momentum 
J = L + 8 {dJ/dt = 0). 

It is also easily checked that 8cs ■ L is a conserved quantity under the above, approximate evolution equations. 
Therefore, as emphasized by Damour Q , when only spin-orbit terms are included the orbital dynamics can be reduced 
to a simple "radial Hamiltonian" H{R,Pfi) = H{R, P^j, = const., S'es • L = const.) describing the radial motion. 
Here P^ = Pi is the momentum canonically conjugated to R ({i?, Pr} = 1). In this case there exists a class of 
spherical orbits satisfying 

R f P n dH{R,PR = Q,L\L-8,s) 

i? = const. , Pr = 0, — =0. (2.31) 

F. Characteristics of the last stable spherical orbit (LSSO) 

When spin-spin interactions are included, those spherical orbits no longer exist as exact solutions. However, as 
spin-spin effects are always smaller than spin-orbit ones, one expects that the above spherical orbits will play the same 
important role as the usual circular orbits in the non-spinning case. In particular, the Last Stable Spherical Orbit 
(LSSO) should play the important role of delineating the transition between adiabatic inspiral and plunge. 

The LSSO for the spinning conservative dynamics is determined by setting 

^-0-^ (2 32) 

where i?o(i?, Pr, P^,---)= H{R, 0, P^, • • • ). 

The physical characteristics of the LSSO for (aligned or anti-aligned) spinning configurations were studied in detail 
m Ref. M within the more fully resummed Kerr-like EOB-Hamiltonian introduced in that reference. It was also shown 
in that the predictions from the latter Kerr-like EOB-Hamiltonian were in good agreement with the numerical 
results on corotating black hole (BH) binaries obtained by means of the helical Killing vector approach [2J|. See 
Table I of . [Note that the agreement with EOB is better when one considers the 3PN accuracy.] The latter Table 
also shows that the numerical results on irrotational BH binaries obtained by means of other approaches based on 
considering only the initial value problem (e.g. significantly differ both from the numerical helical-Killing- vector 
(HKV) results, and the analytical EOB-ones. For recent work improving the numerical implementation of the HKV 



7 



approach (which is closely related to the "conformal thin-sandwich" decomposition) and confirming that it yields 
results that agree well with the EOB approach, see 112, • As all the currently published numerical estimates of the 
physical characteristics of close binaries made of spinning BH's (such as use initial- value-problem approaches 
rather than the physically better motivated HKV one, we shall not try to compare here the analytical EOB predictions 
for spinning configurations with numerical results. On the other hand, it is interesting to compare several different, 
PN-rooted, analytical approaches in their predictions for the binding energy of close BH binaries. 

The most straightforward PN-based analytical approach to the physical characteristics of close BH binaries consists 
of starting from the fully PN-expanded Hamiltonian (|2.15l) . considered as defining an exact dynamics, and then to 
deduce from it the energy and angular frequency of spherical orbits. [We consider here only configurations where 
the spins are parallel (or antiparallel) to the orbital angular momentum, so that it makes sense, even in presence of 
spin-spin interactions, to consider spherical (and actually circular) orbits.] The binding energy of such "PN-expanded" 
spherical orbits is plotted in the left panel of Fig. ^ as a function of the orbital angular frequency, for equal-mass BH 
binaries. As we see from this Figure, the straightforward PN-expanded Hamiltonian does not exhibit any minimum 
of the binding energy, i.e. does not lead to any Last Stable Spherical Orbit (LSSO) in the non-spinning, or aligned 
maximally spinning cases. As the best current numerical results on BH binaries clearly indicate the existence of such 
LSSO's, this disqualifies the use of the fully PN-expanded Hamiltonian (|2.15|) for describing close binaries. 

It has, however, been pointed out that more reasonable results, close to the numerical HKV results, can 
be obtained by plotting, instead of the prediction coming straightforwardly from 12.15|l . the PN-expansion of the 
analytically computed function E{Q) giving the binding energy E as a function of the orbital frequency Q. Indeed, 
one can derive from H2.15II the following explicit PN-expansion of the (invariant) function E{Q) Hj, ji^: 
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(2.33) 
(2.34) 



These functions are plotted in the right panel of Fig.^ Visibly, they are much better behaved than the results plotted 
on the left panel, which came directly from the PN-expanded Hamiltonian. They are also close to the numerical HKV 
results 01- The fact that two expressions, which can both be called "PN-expanded", and which are supposedly 
equivalent "modulo higher PN terms", lead to physically markedly different results lead to conclude that the PN- 
expanded Hamiltonian cannot be used to describe the transition from adiabatic inspiral to plunge. Let us also recall 
that if we consider, in absence of spins, the test-mass limit 77 — > 0, instead of the equal-mass one 77 — > 1/4, the 
PN-expansions (|2.33|) . H2.34|I have been shown (see to be quite inaccurate representations of the known exact 
expression of the function E(yi). Indeed, the 2PN-accurate function (|2.33|l predicts in this limit an LSO frequency 
which is 82% larger than the exact one, while the 3PN-accurate one (|2.34f) predicts an LSO frequency 27% larger 
than the exact one. 

By contrast with these problematic features of PN-expanded results ^, the EOB-approach leads to uniformly better 
behaved results (even if we use it not in the Kerr- like form advocated in Q , but in the form H2.16|l used in the present 
paper). We show in Fig. [21 the EOB analog of Fig. ^ i.e. the function E{Vl) deduced from the EOB Hamiltonian 
(|2.16l) in the (anti-)aligned case. In this EOB case, we have none of the problems entailed by the "PN-expanded" 
approach, and the uniquely defined curve E{Vt) was shown in to agree well with the HKV numerically determined 
curve (for corotating holes). Note, however, that aligned maximally rotating holes lead to a curve which, especially 
in the 3PN case, reaches a minimum (not shown on Figure [2J) only for a rather high angular velocity. 

This property of the aligned configurations (as well as the significant difference between the 2PN-E0B result and 
the 3PN-E0B one) was already emphasized in As it will be important for the present paper, we study it further 
by plotting in Fig. O the dependence on the L- projected effective spin parameter 
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(2.35) 



^ Let us recall here that, in order to be able to describe the transition between inspiral and plunge, we cannot use just the function E{Q), 
but we need a full description of the binary dynamics. Therefore, if we wanted to confine ourselves to a "PN-expanded" approach, we 
would have to use either the PN-expanded Hamiltonian 12.151 . or the corresponding (appropriately truncated) PN-expanded equations 
of moti on. The left panel of Figure 111 shows that this would not be a reliable thing to do. This is also clear from some of the results of 
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(where Scs was defined in Eq. H2.17|l above) of the binding energy E and the angular frequency at the LSSO, i.e. at 
the minimum of the E{Q) curve. This Figure shows four results obtained for equal-mass and equal-spin configurations 
within the EOB approach: (i) the result obtained from the Hamiltonian (|2.16|l when using the 2PN-accurate orbital 
EOB Hamiltonian, (ii) the result obtained from (|2.16l) when using the 3PN-accurate orbital EOB Hamiltonian, (iii) 
the result obtained from the 2PN-accurate Kerr- like Hamiltonian introduced in and (iv) the result obtained from 
the 3PN-accurate Kerr-like Hamiltonian introduced in . [The latter two Hamiltonians are referred to in the caption 
as "SO Resummed" , because they include a resummation a la EOB of the spin-orbit interactions.] In addition, as we 
cannot show on this plot the minimum of the E{Q) curve deduced from the PN-expanded Hamiltonian, because the 
left panel of Fig. 13 shows that it does not exist, we show instead, for comparison purposes, the characteristics of the 
minimum of the PN-expanded functions (|2.33l) . H2.34|l [i.e. the right panel of Fig. (QJ]. 

It is interesting to note on Fig. O that the effect of resumming (a la Kerr) or not the spin-orbit interaction seems 
to be rather small. We see also that, when considering anti-aligned configurations, all calculations give very similar 
results. This is not surprising as the attractive {Hso < 0) nature of the anti-aligned spin-orbit (and spin-spin) 
interaction has the effect of pushing the LSSO outwards, i.e. toward larger-radius, lower-frequency, less bound and 
therefore less relativistic configurations. On the other hand, working at the 2PN or the 3PN level induces, as already 
pointed out in i, a very significant difference for the LSSO characteristics in the aUgr^ed case (positive xl)- In 
this case, because of the repulsive(Hso > 0) nature of the aligned spin-orbit (and spin-spin) interaction, the LSSO 
is drawn towards closer, higher-frequency, more bound and more relativistic configurations. For such very compact 
configurations the repulsive sign (a4 > 0) of the 3PN contribution to the effective potential A{q) further amplifies, 
by a "snow-ball effect" , the tendency toward closer, and more bound configurations. We think that this could be a 
physically real effect due (as confirmed independently by Refs. 0,0]) to the large positive value of the crucial 3PN 
coefficient entering 04, Eq. (|2.14|) . This large positive value for 04 is also needed to improve (with respect to the 2PN 
case) the agreement between the HKV corotating results and the 3PN EOB ones '2^ . It would be interesting to have 
numerical HKV studies of the LSSO for moderately- and fast-spinning aligned BH's to test the predictions made by 
the EOB approach. [The less reliable numerical results of the initial- value-problem of Ref. jij] , which extend up to 
mildly positive values of xl ^ 0.17 are in rough qualitative agreement (especially for the dependence of Qlsso on 
Xl) with the EOB predictions, but their quantitative agreement is too poor to reach a firm conclusion).] 

Let us note in passing that the significant dependence of the LSSO frequency on the effective spin parameter xl 
makes it desirable for data-analysis purposes, when one is content with using adiabatic templates IS], to use at least 
templates whose ending frequency is not fixed say to the usually considered Schwarzschild LSO, but varies with masses 
and spins as suggested by the EOB approach. 

Finally, another consequence of the significant dependence of the LSSO frequency on the effective spin parameter 
XL, is drawn in Fig.0] where we compare the signal-to-noise ratios (SNRs) as function of the binary total mass for 
an optimally oriented equal-mass binary at 100 Mpc. We use LIGO design sensitivity noise curve The SNRs 

are obtained observing the inspiral from 40 Hz until the LSSO predicted by the EOB approach at 3PN order. The 
three curves refer to non-spinning binaries and binaries with xl = —0.875,0.25. Figure ^ reveals a bias towards first 
detecting mostly aligned spinning binaries with high masses, as pointed out in 



III. RADIATION REACTION, INCLUDING SPIN-EFFECTS 



The previous section has reviewed various ways of describing the conservative dynamics of binary systems (including 
spin effects). In the present section, we discuss the inclusion of radiation reaction effects, with emphasis on determining 
the spin-modifications of radiation reaction. Within the Hamiltonian approach, that we use here, radiation reaction 
can be incorporated by modifying the usual Hamilton equations in the following way 

dX' r , . dH 



Here, Fi denotes a "non-conservative force" , which is added to the evolution equation of the (relative) momentum 
to take into account radiation-reaction (RR) effects. This Radiation-Reaction (RR) force F depends, a priori, both 
on the (relative) orbital variables X, P and on the spin variables Si, 82- In the present paper, our aim will be 
limited to determining Fi under the following two simplifying assumptions: (i) we consider only quasi-circular orbits, 
and (ii) we shall retain only the leading spin- dependent terms. After the completion of the work reported in this 
section, there appeared a work of Will 34| dealing with spin-dependent radiation reaction effects in general orbits. 
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As the derivations are not the same, and yield results which we have checked to be equivalent (for circular orbits), 
but expressed in different variables (Hamiltonian X, P here vs. Lagrangian X , V for |3J|), we feel it worth to briefly 
report our derivation. 

Because of the assumption (ii) we look for terms in Fi which are linear in the spin-tensors Sf^ = £ijk [a — \,2) of 
the two considered compact bodies. [Note that spin effects enter the metric only through the spin tensors Sf^, rather 

than through the axial spin vectors S^] As the time derivative of Sfj contains a "small" post-Newtonian factor Gjf? , 
the leading spin-dependent terms in Fi will contain only the undifferentiated spin tensors. 

Using Euclidean invariance, the spin-dependent terms in the force Fi{X,P,Sa) must be a combination of three 
types of contributions: ci Sij , C2 Sij and Ci Sjk X^ P'^, where ci{X, P), C2{X , P) are some scalar functions of 
X, P, while Ci is a vector function of X, P. [Here, Sij denotes one of the two spin vectors. We shall sum over the 
two possible spins at the end.] Imposing that the radiation reaction force Fi be odd under time reversal^ i.e. odd 
under the simultaneous changes X"^ — > X*, Pi ^Pi, Sij ^Sij, tells us that: ci{X,P) must be an even function 
of P, C2{X, P) must be an odd function of P, and the vector Ci{X, P) must be an odd function of P. [Note also that 
Ci must be a true vector, not an axial vector. By parity invariance no Sijk can appear, except in combination with 
S^] Therefore, if we further decompose Ci = C3(X, P) Pi + Ci{X , P) Xi, the coefficient C3(X, P) must be even in P, 
while the coefficient C4{X, P) must be odd in P. 

At this point, our simplifying assumption (i) above (quasi-circular motion) will bring a drastic simplification. Indeed, 
a time-odd scalar must contain an odd power of the combination X^P^. However, this combination vanishes along 
circular orbits (and is therefore subleadingly small along adiabatically inspiralling orbits). To leading order the scalar 
coefficients C2 and C4 therefore vanish, and we conclude that Fi contains only two independent spin contributions: 
ci{X , P) Sij X^ + C3{X , P) Pi Sjk X^ P'' . It will be convenient in the following to further decompose the vector 
Sij X^ entering the first contribution (which is orthogonal to X^) into its component along the direction of Pi, and 
its component orthogonal to Pi, say 

{S,, = {S,k - P^ Pk/P^) Skj X^ = S,, X^ + ^ Sjk X^P" . (3.2) 

It is easily checked that, along circular orbits (AT* Pi = 0), the vector (|3.2() is orthogonal both to P and to X. Therefore, 
{Sij X^)-^ is parallel to the orbital angular momentum (axial) vector 



= Eyfe X^ Pk . (3.3) 

It is easily checked that 

{S,, = ^{p.s)L, = ^{P-S)L,. (3.4) 

Finally, adding the usual spin-independent radiation reaction (parallel to Pi for circular orbits), and summing over 
the two bodies, we conclude that the RR force can be written as 

F,{X,P,Sa) = BP, + Aa{St,X^)^ =BP + ^[P-Sa)L,, (3.5) 

a=l,2 a=l,2 

with 

B = Bo + ^ Ca A^P^ = Bo + ^ Ca i • 5a , (3.6) 

a=l,2 a=l,2 

where Bo, Ca and Aa are some time-even functions of -X" and P. 

To determine the coefficients Bq, Ca and Aa, we now impose that there be a balance between the losses of mechanical 
energy and angular momentum of the system due to the additional force Fi in the Hamilton equations of motion H3.1() 
and the losses of energy and angular momentum at infinity due to the emission of gravitational radiation. Let us first 
recall that, in the Hamiltonian formalism, the quantities 

E{t)=H{X{t),P{t),Sait)), 



Jy (t) = X' P, - X^ P, + si + si , 



(3.7) 
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are exact constants of the motion in absence of RR force in Eqs. H3.1() and in the Hamiltonian equations for spin 
evolution. When adding the RR force F in Eqs. (|3.1|) (and no corresponding RR torque in the spin evolution 
equations) we find that E and J evolve as 

dE_dH ■ 

^ ^ X' F, - X' F, . (3.9) 

Inserting in Eqs. (|3.8I3.9|I the expression H3.5fl for the RR force F^, we can easily evaluate the averaged losses of energy 
and angular momentum. Along (quasi) circular orbits the various scalar coefficients -Bq, C'a, Aa are time-independent 
(because all basic scalars, X^,P^,X • P = 0, are constant). One then finds that dE/dt is time-independent, while 
dJij/dt depends on the orbital phase only in spin-dependent terms and through the tensor X' Pj. Decomposing the 
latter tensor into 

X^ P, = i(X' P, + X^ P,) + ^{X^ P, X^ P,) A Q X-'-) + i , (3.10) 

one easily sees that its orbital average is simply {X"^ Pj) ~ [We consider averages over the orbital period, 

considering all more slowly evolving quantities, such cis fixed during one orbital period.] 

When evaluating Eq. (|3.8() along circular orbits, we cannot use the Newtonian approximation X^ ~ P'/m because 
we wish to obtain the coefficient B with a high post-Newtonian accuracy. However we can instead use JC^ Pi — (j> P^ — 
uj\L\ where oj — (f) — V/ R denotes the orbital angular frequency. We finally obtain 

dE 

— =Bu;\L\, (3.11) 



dJ 



dt 2 

a=l,2 



where A = L/\L\ denotes the unit vector along the orbital angular momentum. Note that Eqs. H3.11|l . H3.12|l predict 
a link between energy loss and angular momentum loss, namely 

dE /dJ\ , ^ 

-J7-^>^-(-J7) ■ (3-13) 



dt \ dt 



To obtain the values of the coefficients B and Aa (and to test the prediction (13.131) ). we need to compare Eqs. H3.11|l . 
(|3.12l) with the values of the averaged fluxes of energy and angular momentum at infinity. The spin contributions to 
the latter losses have been computed by Kidder . However, one must be careful with the fact that Kidder expressed 
most of his results in terms of harmonic coordinates, with the choice of a covariant spin supplementary condition: 

First, using the results of Ref. 9] as they are, one straightforwardly checks that the relation H3.13|l is satisfied. This 
is a check that it is enough to include RR effects in the orbital equations of motion (|3.1|l , without modifying the spin 
equations of motion. [For a direct dynamical check, see [s^.] Then, to obtain the value of the coefficient B we can 
simply use the result 13.11(1 . namely 

u} \L\ dt 

where it remains to express dE/dt (computed as a fiux at infinity, using Ref. (9|) in terms of our basic (Hamiltonian) 
dynamical variables. One way to proceed would be to transform the harmonic-coordinates results of into ADM 
coordinates (with the corresponding spin condition Siq -I- \SijV^ = [30l|'). The transformation linking the two 
coordinates has been worked out in j30j (for the spin-dependent terms) and in ^] for the spin-independent parts. 
However, a simpler way to proceed is to eliminate references to specific coordinates by expressing dE/dt (for circular 
orbits) in terms of the gauge-invariant orbital frequency uj. Adding also, for better accuracy the recently completed 
3PN flux contribution 50, 51, 52], we have 

^ = -y + hiri) vl + [/3(r/) + /aso] vl + [^(r]) -h Ass] + f,{rj) vl + + f,, ln(4«^) -I- /^(r;) vl} , 

(3.15) 
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where = (GMw)-"^/^, where the spin-independent flux coefRcients /2(??), • ■ ■ , f7{v)i *re given by 

1247 35 



336 12 



(3.16) 



(3.17) 



fM = - 



44711 9271 



9072 



65 



W^+18^ 



/5(^) = 



/8191 583 
= - -r^ + —1 



V 672 



24 



77 TT, 



(3.18) 



hiv) 



6643739519 16 



69854400 



1712 
T05 



7b 



134543 
7776 



41 
48 



94403 
3024 



775 

324 



(3.19) 



1712 
T05" 



(3.20) 



f7{v) = 



16285 214745 



504 



1728 



193385 
3024 ' 



with jE being Euler's gamma, and where the spin-dependent corrections to the latter flux coefficients are 





5 






(t- 


5 


mi ^ 


A • S2 




^4 


mi J 


GM^ 




^4 


m2 J 


GPP 



(3.21) 



(3.22) 



Ass 



[289(A ■ Si)iX ■ S2) - 103 Si ■ S2] + 0(31) + 0(31) 



(3.23) 



The present work was aimed at determining the leading spin-dependent terms, i.e. the ones linear in 5'i and ^2, as 
exemplified in the correction fsso, Eq. (|3.22|) . to the coefhcient /s — An. For completeness, as the link H3.14|l between 



the "longitudinal" part of RR, 



long 



B Pi, and the energy loss, is clearly general, we have also used Kidder's results 



to write down the part of B which depends on the product S\ S'^ . The numerically similar contributions which 
depend on S\ S{ and ^2 5'^ have not yet been determined. Only partial results are known. For instance, Poisson [ss" 
has derived a contribution to /4gs of the form 



3{X -31)^-31 , 3{\-S2) 



G^m^M^ 



(3.24) 



but many other additional contributions 0{Si) + 0(82) have not yet been computed. 

Let us finally turn to the determination of the other spin-related coefRcients, Aa, in Eq. H3.12|l . Again, we have 
the technical problem that Ref. expressed its results in terms of harmonic-coordinate quantities. Namely, Eq. 
(4.11) of Ref. PI expresses the total angular momentum loss dJ/dt in terms of the harmonic distance r and of the 
harmonic-coordinate "Newtonian orbital momentum" Ln = ^.x x v (where x and v denote the relative harmonic 
position and velocity). A simple way to convert this result to our ADM distance R and our ADM total orbital 
momentum L = X x P is to relate Lj^ to L hy comparing the expression (4.7) of Ref. for the (gauge-invariant) 
conserved total angular momentum J with the corresponding simple ADM expression H3.7|l . This yields a result of 
the form 



■IN 



cL 



a=l,2 



Xa Sa 



(3.25) 



where the coefficient c is not needed for our present purpose, and where, following the notation of 9], Sa = Xa ^a- 
Inserting Eq. (|3.25|) in Eq. (4.11) of Q allows one to compute easily the part of dJ/dt which is proportional to the 
projection of Sa orthogonally to L. 
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This yields the following expression for the coefficients Aa in Eq. H3.5|l 

where a' ^ a (e.g. a' = 2 when a = 1). Summarizing, the radiation reaction force to be added to the Hamiltonian 
equations of motion (|3.1|) reads 

where the energy loss (expressed in terms of the orbital frequency lu, or equivalently of Vi^ = {GMcoy^^, and of the 
spin variables) is given by Eqs. H3.15|l - H3.23|l . We have checked that, after taking into account the relation between 
the Hamiltonian variables X, P and the Lagrangian ones X, V (which involves spin-dependent terms because of the 
first Eq. H2.23|l . Eq. H3.27|) agrees with the circular limit of Eq. (1.6) of (3^] (which assumes the same spin condition 
as we do) . 

Refs. 01 USES have shown that (at least in the test-mass limit where one can compare analytical and numerical 
estimates) it is generally advantageous to replace the Taylor series in curly brackets on the right hand side (R.H.S.) of 
Eq. H3.15|l by its (suitably defined) Pade resummation. In particular. Porter and Sathyaprakash have compared 
"Taylor" and "Pade" approximants for the flux function of a test particle around a Kerr black hole with the exact 
numerical estimates [s^ l and concluded that Pade approximants are, when considering all values of the spin parameter, 
both more effectual (i.e., larger overlaps with the exact signal) and more faithful (i.e., smaller biases in parameter 
estimates) than Taylor approximants. [We use here the terminology introduced in 19].] In view of this, and as 
was already advocated in Q, we consider that that the best way to incorporate a radiation reaction force in the 
EOB approach is to insert Pade approximants of the flux function (R.H.S. of (|3.15l) ) in (|3.27|l . However, for added 
generality, we shall also consider the case where we leave the flux function as a plain Taylor series. Note that, when 
considering arbitrary values of the dimensionless spin parameters for the two holes xi = Si/Gm\, xi = SilGm\ we 
used the normal "direct" (i.e., lower-diagonal) Pade-approximants, instead of the "inverse" (i.e., upper-diagonal) ones 
used in 55]. For some values of the spin parameters both the lower and upper diagonal Pade-approximants have poles. 
When this occurs, we apply the Pade-approximant only to the non-spinning part of the flux and add the spinning 
terms separately. There exist other Pade-approximants in which poles are absent and it would be very desirable to 
determine them in the entire parameter space. This is beyond the scope of this paper. 

In Figs.|31|Hlwe show T- and (lower-diagonal) P-approximants at 3.5PN order for an equal-mass binary and several 
values of the dimensionless spin parameter x = Xi = X2- We notice that the T- and P-approximants are much closer 
in the anti-aligned cases than in the aligned one. Since the calculation of the non-spinning flux function at 3.5PN 
order has been completed only recently M^, in Fig.|7|we contrast the T- and (lower-diagonal) P-approximants at 3PN 
and 3.5PN order for an equal-mass non-spinning binary. 



IV. DEFINITIONS OF THE INITIAL AND ENDING CONDITIONS OF TWO-BODY MODELS 

As clear from the comparison of the left and right panel of Fig. because of the bad behaviour of the PN 
Hamiltonian near LSSO, we propose as our best bet, for describing in a physically reliable manner the non-adiabatic 
evolution of BH binaries, and their transition between inspiral and plunge, to use an EOB-resummed Hamiltonian, 
we shall, for more generality, consider, and compare, in this Section several types of two-body models. 

To define a specific model we must make various choices: (i) choice of a PN-expanded (or "Taylor-expanded") 
Hamiltonian (say "TH") versus an EOB-resummed Hamiltonian (say "EH"); (ii) choice of a Taylor-expanded flux 
function (say "TF" ) versus a Pade-resummed one (say "PF" ), and finally, (iii) choice of the PN accuracies used both in 
the Hamiltonian (say n PN) and the flux function (say to PN). This leads to models denoted, for instance, THTF(n, to), 
EHTF(n, to), EHPF(n, to). In addition, as we are here mainly considering the evolution of spinning binaries, we shall 
add an initial letter S to recall that fact. This leads to models denoted as STHTF(n, m),..., SEHPF(n, m). To 
simplify, we shall only consider the PN accuracies (2, 2.5) or (3, 3.5). To further simplify, we shall focus on comparing 
"fully Taylor" models (i.e., STHTF), to "fully resummcd " ones (i.e., SEHPF). Finally, this leads us to considering 
only four models: STHTF(2, 2.5), STHTF(3, 3.5), SEHPF(2, 2.5), and SEHPF(3, 3.5). [Note, as discussed above, that 
because of the appearance of spurious poles in a few tests, we applied in those cases, the Pade resummation only to 
the non-spinning part of the flux.] 

An important parameter in our present model building is to decide when to stop the evolution. This issue was 
already tackled in Ref. 7]. There, because we were using an EOB Hamiltonian, and were considering non-spinning 
BH's, we found that we could follow the "plunge" (after LSO crossing) down to a (Schwarzschild-like) radius ~ 3M, 
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at which point we could match to a ring-down signal made of least-damped quasi-normal modes. It was found in [7| 
that, contrary to what the usually employed word "plunge" suggests, the inspiral motion after crossing the LSO was 
staying "quasi-circular" , with a kinetic energy in the radial motion staying small in absolute value, and smaller than 
0.3 times the kinetic energy in the azimuthal motion down to i? ~ 3M. In our "spinning" evolutions the situation 
is more complicated (notably when considering large aligned spins, and also, for evident reasons, when considering 
Taylor-expanded Hamiltonians). We leave to future work a detailed discussion of the matching to ring-down. We 
decided to stop the evolution as soon as one of the following inequalities ceased to be fulfilled: 



\R\ < 0.3|Vt|, (4.1a) 
P^/B(i?) < 0.3Pl/R\ (4.1b) 

\EnR\ > 0.1|i^Sri> (4.1c) 
R > aM , (4. Id) 

where B{R) = D{R)/A{R) [see Eqs. (glSJ, ^T^]. Criteria (|Oa|) - l|4TH|l ensure that the evolution does not extend 
too much beyond circularity, on which our formulation for radiation-reaction force relies. The quantity Vt is the 
tangential velocity (i.e., orthogonal to the relative separation vector X). Criterion (|4.1c|l is used to avoid going into 
regimes where the GW energy flux goes to zero (e.g., for Taylor-expanded flux at 2.5PN order). Criterion Ij4.1dp (in 
which a ~ 1 when using the ADM-coordinate Taylor-expanded Hamiltonian, and a ~ 2 when using the Schwarzschild- 
like EOB Hamiltonian) terminates the evolution at a very small radius, in case all of the above criteria fail to take 
effect. 

In all cases, the instantaneous GW frequency at the time when the integration is stopped defines the ending 
frequency for these waveforms. We shall also consider below extended waveforms obtained by matching a ring-down 
signal when this ending frequency is reached. 



A. Initial conditions: quasi-spherical orbits 

[In this section we shall use natural units c = 1 = G and set M — 1] 

In absence of radiation reaction (RR) , spherical orbits with constant radius and orbital frequency exists under spin- 
orbit interactions, but cease to exist when spin-spin interactions are present (except in special situations when the spins 
and the orbital angular momentum are all aligned/anti- aligned). When radiation reaction is treated adiabatically, an 
initial spherical orbit will evolve into a sequence of spherical orbits, due to Eq. (|3.13|) . In this section, we formulate 
a prescription to construct initial conditions for non-adiabatic evolutions, which lead to quasi-spherical orbits under 
spin-orbit interaction. 

With spin terms kept only up to the spin-orbit order, the Hamiltonian can be re-written into a simpler form. 



i/(i?, Pi?, L, Xl) = ffnospin(i?, Pi?, i) + 2 

Here i?nospin a-re terms in the Hamiltonian that do not involve spins, and 



i?3 



(4.2) 



L = |L| 



XL = S 



off 



(4.3) 



In this form, the Hamiltonian depends on four quantites, {P, Pr,L,xl}: in which L and xl both depend on 
{6', (/), Pe, P^}, while xl also depends on the spins. In absence of radiation reaction, the conditions for spherical 
orbits written in terms of partial derivatives (indicated by a subscript i) with respect to the four independent vari- 
ables {P, Pr, L, Xl}, read 



[R\ = Q 









XdRJe. 








= 0. 



(4.4) 
(4.5) 



[Here the subscript c indicates canonical partial derivatives. In the rest of this section, we shall continue to use i and 
c to distinguish between these two types of partial derivatives.] With L and xl being conserved quantities, conditions 
(|4.4|) and (|4.5|) will remain satisfied if they are initially satisfied — which proves the existence of spherical orbits. 

We now construct initial conditions for spherical orbits, in absence of radiation reaction, based on Eqs. H4.4|) and 
(|4.5|) . In numerical evolutions, given a source coordinate frame, {ex,ey,ez}, we specify spherical orbits with the 
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following initial kinetic parameters: the orbital frequency luq: the orbital orientation (i.e., the normal direction to the 
orbital plane [Ln]o = x ^)/\^ ^ ^1): the spins [Si^2]o, and the direction of initial orbital separation N = 
which can in turn be given by an initial orbital phase i^orb, calculated with respect to the reference direction of 
[Stot X Ln]o, 



No 



[Stot X Ln]o 



[Ln]o X [Stot X Ln]o . 
b + - — . sm ( 



|[Stot xL 



[[Stot X Ln]o| 

We calculate initial values for {X, P} in three steps: 

1. We first apply a rotation TZ such that [Ln]o ^ and Nq ex- 

2. In spherical polar coordinates, the above step implies 

(j)o = 0, eo = TT/2. 



(4.6) 



NjO 



(4.7) 



[The (j)Q here should not to be confused with the orbital phase (poib above.] Then, we specify the initial frequency 
ujQ and impose 



UJQ = (j)o 



OH 







dH 
'dR 



= 0, 



(4.8) 
(4.9) 



and solve for the four variables {R, Pr, Pg, P^jo. 

3. Finally, we apply the inverse rotation TiT^ to the entire system, obtaining a set of initial spherical-orbit conditions 
consistent with the specified initial kinetic parameters. 

When radiation reaction is included, we proceed as in Ref. and modify Eq. 1)4. 4() to include a non-zero R, 
according to the prediction from adiabatic evolution. 



{dE/dR) 



sph 



(4.10) 



in order to prevent radial oscillations. [The subscript sph in Eq. (|4.1U|) and below denote quantities evaluated along 
spherical orbits.] Equations (|4.5I) can be kept unchanged, since Pr is second order in radiation reaction. We now 
calculate {dE / dR)sph in terms of the simphfied set of independent variables, {R.Pr.L.xl}- Consider neighboring 
spherical orbits in an adiabatic sequence, we have 



dH = 



dH 
'dR 



dR 



dH 
dP~B, 



dPn 



dH 
'dL 



dL 



f dH 
\dXL 



dXL 



(4.11) 



dH 
'dR 



-0, =0, Pfl -O, 



It is straightforward to deduce that 



dE 
'dR 



( dH 
\~dL 



fd^H\ 



sph 



/ d^H \ 
KdRdLj, 



dH 
dxL 



( dH 
\dL 



dPn = . 



/ d'^H 
\dRdxL 



(4.12) 



/ d'^H \ 
[dRdLj^ 



dXL 
dR 



(4.13) 



sph 



will be ignored 



The second term on the right-hand side can be ignored, as we argue later in this section, because xl is still conserved 
to a high accuracy even in presence of radiation reaction. In special configurations with L = (or equivalently 
9 = 7r/2, Pg = 0) we can re-write Eq. (|4.13(l in terms of canonical variables in spherical-polar coordinates: 



dE 
'dR 



sph 



m/dP^)c{d^H/dR^)c 



(d^H/dRdP^), 



(4.14) 



9=TT/2,Pe=0 
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We also note that when L is known to be e^, we can calculate E'rr, from Eq. (|3.15() right away using only ujq and 

[Sl,2]o- 

We can now construct quasi-spherical initial conditions when radiation reaction is present. As done in Ref. |2f , up 
to leading order in radiation reaction, we only need to augment our no-radiation-reaction initial conditions with a 
non-zero Pr, with initial values for all other canonical variables unchanged. In order to do so, we insert three more 
steps between steps 2 and 3 above: 

2a During step 2, we have obtained a set of spherical-polar-coordinate initial conditions, for a rotated system with 
{No, [Ln]o} — {^XT^z}- The canonical angular momentum, [L]o of this system, though, will not in general be 
along e^. However, being orthogonal to Nq = e^^, it must be within the ey — plane. We now apply a further 
rotation TZ' around Nq = to the entire system, such that afterwards we have {Nq, [L]o} = {ex,ez}, i.e., 
00 = tt/2 and [Pg]o - 0. 

2b Now that Eq. (|4.14() is applicable and Ebr is readily obtainable, we insert them into Eq. (|4.10(l and obtain [R]o. 
[Note that in this process we use the set of initial conditions already obtained for a spherical orbit in absence of 
radiation reaction, with Pr = 0.] From this [R]o , we obtain the initial value of [Pr]o to insert into our existing 
set of initial conditions: 



{R]o 




fdH\l 






KdPnJJ 


Pr^O 



(4.15) 



2c Gathering our new set of {R, 0, 0, Pr, Pg,P^}o, we obtain the Cartesian-coordinate variables, and apply a inverse 
rotation {TZ')^^ to the entire system. [Now again we have {Nq, [Ln]o} = {^x^^z}, and are ready to proceed to 
step 3.] 

A straightforward analysis of the various error terms allows us to conculde that the fractional error of assuming 
that XL = Scff ■ L is constant along the adiabatic evolution is of 3PN order. 

We note that our steps 1, 2, (2a, 2b, 2c), and 3 can still be applied to give initial conditions, even if the Hamiltonian 
contains spin-spin terms, although the orbits that follow will in general have oscillatory radius and orbital frequency, 
due to the non-existence of quasi-spherical orbits. In Fig. [S] we show the evolutions of r/(ruj) with (dark curves) and 
without (light curves) spin-spin terms, for (10 -|- 1Q)Mq (left panel) and (15 -I- 5)Mq (right panel) binaries. We start 
evolution at 40 Hz, with {0s^As^Os2, 4>S2) = (60°, 90°; 60°, 0°), and show the evolution up to 200 Hz. 

V. COMPARISON OF WAVEFORMS AND EVALUATION OF OVERLAPS 

In harmonic coordinates, the gravitational wave emitted by a binary system at the leading quadrupole order, in 
terms of metric perturbation at a distance D, is 

Using the leading-order equation of motion, Xk = —MX^/R^, we re- write the normalized perturbation Hij as: 

H,, - 4^^ {v, V, - M^j^^ . (5.2) 

Here Xi and Vi = Xi can be obtained straightforwardly by solving the Hamilton equations. Depending on the wave 
propagation direction and the orientation of the detector, the metric perturbation hij has to be contracted with an 
appropriate "detection tensor" to give the actually detected waveform. For this we refer, for instance, to Sec. IIIC of 
Ref. [13 and Sec. II of Ref. 15] (in particular see Eq.(15)). 

Following Ref. (IS] , the parameters in precessing binaries can be distinguished in binary local parameters 
{m\,m2, Si, S2,0si,Qs2,4>si — 4'S2\, binary directional parameters {6'l,0l,0si + 4'S2\ (which determine the orien- 
tation of the binary as a whole in space) and directional parameters {Q,ip,0,(j3,ip}, describing the GW-propagation 
and the detector orientation [see Table I in Ref. and discussion around it]. To these parameters we need to add 
the initial time and the initial orbital phase. 

In the precessing convention introduced in Ref. [l^] , the GW signal can be neatly written in terms of: (i) parameters 
depending on the observer's location and orientation, {Q,ip,0,(f),^}, which are time independent, initial time and 
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initial orbital phase (henceforth denoted as extrinsic parameters) and (ii) parameters depending on the details of the 
dynamics, {mi, TO2, Si, S'2, 6'si, ^32, 4'Si ~ 4'S2} (henceforth denoted as intrinsic parameters). The GW signal does not 
depend on the binary directional parameters, {^l, 0l, 0si + 0S2}, since those parameters can be re-absorbed in the 
definition of the source frame at initial time and in the directional parameters {Q, ip, 9, 0, ip} through a rigid rotation 
of the detector-binary system. 

The distinction between extrinsic and intrinsic parameters is due originally to Sathyaprakash and Owen [H^ l57| | . 
Extrinsic parameters are parameters which change the signal shape in such a way that we do not actually need to 
lay down templates in the bank along those parameter directions, saving computational time. By contrast we need 
to lay down templates along the directions of the intrinsic parameters. In Refs. 15, 16], a semi- analytical method to 
maximize over the extrinsic parameters in precessing binaries has been proposed. 

In view of the bad performances of the Taylor-expanded Hamiltonian'^, we a priori expect that the waveforms 
computed from STHTF models will be significantly different from the SEHPF ones. It remains, however, interesting 
to measure their difference in the data-analysis sense, i.e. to compute the overlaps between the two types of waveforms. 
If it happened that, after maximization over all possible parameters, the overlap between the two types of signals 
were very close to unity, one could still consider the Taylor models as effectual (in the sense of T^) representations 
of the EOB models. However, for practical reasons, we did not try to embark on a full maximization of the overlaps. 
For simplicity, we only tackled the maximization over the extrinsic parameters, and not on the intrinsic ones. The 
resulting partially maximized overlap is therefore only a lower bound of the fully maximized overlap. Still, this result 
can be considered as a reasonable measure of the "closeness" between the two sorts of models (especially because we 
do not wish to use models which would have significantly different physical parameters). 



A. Lack of "closeness" between Taylor and Effective-One-Body models 

In Tables III and HTl we study the closeness (in the sense just defined of overlap maximized only over the extrinsinc 
parameters'') between STHTF(3,3.5) and SEHPF(3,3.5), as well as between STHTF(2,2) ^, and SEHPF(2,2.5) models. 

We consider three typical binary masses (10 -I- 10)Mo, (15 -I- 15)Mq and (15 -I- 5)M0, and several initial spin 
orientations ^ . We always fix the pattern functions — \, — Q and GW propagation parameters = 7r/4 
and 1^ = [for notations and definitions see Sec. IIIC of Ref. 16] and Sec. II of Ref. "l^]. The initial frequency is 
always set to /in = 30 Hz and the ending frequency /end is determined by one of the criteria in Eqs. H4.1a|l - lj4.1d|l . 
In Tables III and ITU the two black holes are assumed to carry maximal and half- maximal spins, respectively. We list 
the ending frequency, the LSSO frequency and the BH radial separation at ifin for the template and target, together 
with two types of overlaps: the overlaps maximized over the initial time and initial orbital phase only (pmax,2), and 
the overlaps maximized over those parameters and {0,(p,a = f{9,(j),ip)}, as well, (pmax.s), using the semi-analytical 
method suggested in Ref. O. Table jlj and jnjalso contains the non-spinning case. 

As these tables show, the two types of models are not at all "close to each other". The overlaps are indeed 
quite low, as low as 0.32. The overlaps evidently increase when we maximize over five rather than two extrinsic 
parameters, but not dramatically, and only for binaries with high and comparable mass, with initial spins lying in 
the half-space opposite (with respect to the orbital plane) to the initial orbital angular momentum. In this case the 
dynamical evolution is shorter, since the LSSO occurs at lower frequency [see also Figs. 13 0], and the differences in 
STHTF(3,3.5) and SEHPF(3,3.5) can be compensated by an offset in the extrinsic parameters of the template with 
respect to the target. Moreover, both the conservative dynamics for circular orbits and the GW flux predicted by 
SEHPF-approximants and STHTF-approximants, are closer in the anti-aligned case than in the aligned case, as can 
be see in Figs. jSj and 

However, when the binary mass ratio is significantly different from 1, the number of modulational cycles increases, 
and the differences in the two models can no longer be compensated by re-adjusting the template extrinsic parameters. 
When the initial spins are lying in the same half-space of the orbital angular momentum, the evolution is longer, the 



^ We have in mind here the absence of LSSO. Recall also that when the binary mass ratio is significantly difi'erent from one, one can firmly 

conclude that the Taylor-expanded Hamiltonian is a poor representation of the dynamics, while the EOB-resummed one is a good one. 

More precisely, we do the maximization over the extrinsic parameters of the EOB model. Though this introduces an asymmetry in the 

definition of the closeness, we do not expect this asymmetry to be physically significant. 
^ We use SHT(2,2) instead of S TH TF(2,2.5) because for equal-mass binaries the Taylor-approximant for the fiux at 2.5PN order becomes 

negative for large values of v llHl 

® For these data we always refer the initial spin directions to the initial direction of the orbital Newtonian angular momentum, as specified 
in Fig. 4 of Ref. [ij, and we set the initial direction of the Newtonian orbital angular momentum along the x-axis of the source frame 
(i.e., wc fix = n/2 and (j>i, = 0, see Fig. 3 in Ref. IT^ I. 
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LSSO happens at high frequency [see also Figs. Ellll], and in this case, even for high, comparable masses the differences 
both in the conservative and non-conservative late dymanics in the two models cannot be compensated by a bias in 
the template extrinsic parameters. 

From Tables Ull we notice that all the above considerations apply also at 2PN order, where the differences in the 
conservative and non-conservative dynamics of STHTF and SEHPF approximants are even larger. We checked that 
these considerations do not change much when spins are smaller, say half-maximal. 

Having confirmed that Taylor models cannot be considered as being effectively close to the FOB ones, we shall only 
use in the following the a priori better FOB models. 

B. Negligible influence of the of the "transverse component" of the radiation reaction force. 

Having in mind possible simplifications of the models, we first investigated the relevance of the second term in the 
R.H.S. of Eq. (|3.27() . i.e. the component of the radiation reaction force which is "transverse", in the sense of being 
directed along L, and therefore orthogonal to the main "longitudinal term", which is parallel to the momentum P). 
In Table IhH we study the influence of this transverse component of the RR force on the dynamics and the waveforms. 
For the binary masses (10 -I- 10)Mq, (15 + 5)Mq and a few initial spin orientations, we evaluate the same quantities 
of Tabled when including and not including the RR force along L [see second and third term in Eq. H8.27|l ]. We give 
here only the overlaps maximized over five extrinsic parameters. We find that Pmax.s is larger than ^ 0.98 in all cases. 
We therefore conclude that it would suffice to use a simplified RR force parallel to the linear momentum P (as in the 
non-spinning circular case). One should, however, include, for better accuracy, in the coefficient of Pi in Eq. (|3.27() 
the spin-dependent terms. 

C. Influence of the resummation of the "longitudinal" part of the radiation reaction 

We consider here the "longitudinal" part of the radiation reaction, i.e. the first term on the R.H.S. of Eq. 13.27|l . 
This term is given by the fiux function, which was written in Eq. (|3.15|l above as a straightforward PN-expansion. 
One can therefore either leave this longitudinal component in non-resummed, "Taylor" form, or choose to resum it 
by means of Pade approximants. In Table ITVl we investigate how the choice of the flux function (Taylor-expanded or 
Pade-resummed) may affects the dynamics and the waveforms. Using in all cases the FOB Hamiltonian to describe 
the dynamics, we evaluate the overlaps between models using a Taylor fiux and models using a Pade fiux. [We 
maximize over the five extrinsic parameters of the FOB model.] We find that, when the initial spins are lying in the 
same half-space of the orbital angular momentum, after maximization over the five extrinsic parameters, the overlaps 
are reasonably large (larger than 0.84), but still lower than unity. We obtain much higher overlaps when the initial 
spins are not lying in the same half-space of the orbital angular momentum. These results are consistent with Figs.EI 
and Fig. 

If we assume that the equal-mass flux function is a smooth deformation of the test-mass limit one, since previous 
findings [3, 0, [s^ in the test-mass limit case pointed out the usefulness of Pade-resumming the fiux function, 
we would conclude that Pade-resummed fluxes are better approximants of the numerically determined flux also in th 
equal-mass case. 

D. Negligible influence of the quadrupole-monopole terms 

Still in the spirit of trying to simplify the models to their crucial elements, Table Ivl investigates how waveforms are 
affected by the quadrupole-monopole terms, and Table IVTl studies how the evolution obtained by averaging the spin 
terms over a period may differ from the non-adiabatic evolution. Considering the high values of Pmax.s we obtain in 
both cases, we can say that the quadrupole-monopole interaction and the adiabaticity of the spin terms, have little 
physical effects over the dynamics and waveforms. The differences can be compensated by re-adjusting the template 
extrinsic parameters. 

E. Influence of the initial orbital phase 

Finally, we investigated the influence of the initial orbital phase (all other quantities being flxed) on the waveform. 
In an adiabatic evolution in which spin terms are averaged over a period the joint evolution of Ln and S is not 
affected by the initial orbital phase. As a consequence, two configurations with the same initial values for Ln and 
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S, but different orbital phases (po-ch will keep the difference between orbital phases unchanged through the evolution. 
This may not be true in our non-adiabatic evolution for two reasons: (i) the spin-spin interaction Hamiltonian 
depends explicitly on the separation vector N [see Eq. H2.18|l ]. and (ii) the evolution depends on the canonical orbital 
angular momentum, which is not orthogonal to the orbital plane. We illustrate this feature by evolving two maximally 
spinning (15 + 15)M0 binaries with initial orbital phases (at 40 Hz) differing by 7r/2, and all other parameters identical: 
{9si,(t>Si',0s2,4'S2) — (60°, 90°; 60°, 0°). In Fig. El we plot the difference A(/)oib between their relative orbital phases 
measured with respect to Ljv- This difference grows in time, and accumulates around 270° by the end of the evolution. 
We also show waveforms detected with {F-^- , Fy ; Q , ip) = (l,0;7r/4, 0) in Fig. EH Their phases start out to differ by 
180° as expected, and non-adiabatic effects drive them away by more than 2 cycles toward the end of the evolution. 
We notice that comparing the waveform is less straightforward than comparing the relative orbital phase, because the 
waveform phase can differ from twice the orbital phase, due to precessions. 



VI. LOSSES OF ENERGY AND ANGULAR-MOMENTUM AND THE WAVEFORM INCLUDING 

RINGDOWN 

In the following, we use as model the "best bet" we can make, i.e. the spinning FOB Hamiltonian^ with a Pade- 
resummed flux. Both being taken to the highest PN-accuracy available, i.e. n — 3,m — 3.5, in the notation used 
above. 

In Ref. 58], using the non-spinning EOB Hamiltonian at 2PN order, it was found that the energy emitted during 
the plunge is ~ 0.7% of M, with a comparable energy loss ~ 0.7% of M during the ring-down phase. This gives a 
total energy released beyond the LSSO in the non-spinning case of ~ 1.4% of M to be contrasted with 4 — 5% of M 
estimated in Ref. |23|, where the authors use a combination of numerical and perturbative approximation methods. 
Note also that Flanagan and Hughes predicted 10%M for inspiral and plunge, and ~ 3%M for ring-down 
phase. 

Here, to have more confidence in our EOB-based estimates, we decided to use three different ways of evaluating 
the energy radiated in the spinning case. We used, at once, (i) the change, along the evolution (between some initial 
frequency and some final one) in the numerical value of the Hamiltonian H2.16|l {6 Eh), (ii) the time-integral of the 
square of the third derivative of the quadrupole moment Itj with i,j = 1, 2, 3, i.e. 



^^i[x,XJ- -s.jX'^ Xk . (6.1) 



dt 5 dt^ dt^ «j p \^ » J 3 y 

and (iii) the time-integral of the energy flux carried away by our leading-order quadrupole waveform, 

with HJ^ being the trace- free part of the normalized metric perturbation, Hij [see Eq. 1)5. 2|l above]. 

In Fig. ^] we compare the accumulated energy release from these three prescriptions, for a (15 -I- 15)M0 maximally 
spinning binary with a generic set of spin orientations when starting evolution at /gw = 30 Hz: [Os^ , (jis-i ; ■, 4>S2) = 
(60°, 90°; 60°, 0°). In the left panel, we keep radiation reaction force at the Newtonian order, while we use 3.5PN 
Pade flux in the right panel. As we see from the figure, these prescriptions differ more from each other when 3.5PN 
radiation reaction is used instead of Newtonian — this is consistent with the fact that both 6Ei and 5Eh involve 
quadrupole radiation only; furthermore, for lower frequencies the 5Eh curve lies below those of 5Ei and 5Eh, which 
is consistent with the fact that Post-Newtonian GW luminosity is in general smaller than the Newtonian prediction. 
The difference among 5Ei and 5Eh can be attributed to the difference between PN (in our case EOB at 3PN) and 
Newtonian dynamics, which seems to be small till around /gw = 200 Hz in our case (which corresponds to v k, 0.45). 

The rather satisfactory agreement between the various ways of estimating the energy loss gives us some confidence in 
our EOB-based estimates. In the following, we shall use the a priori best estimate (because it is the one which involves 
the highest PN accuracy): the one based on the change in the total EOB Hamiltonian H. [We use Pade-resummed 
fluxes, and the two combined PN-accuracies {n,m) = (2,2.5) and {n,m) = (3,3.5).] 



^ We did not investigate the "closeness" between the models derived from the spinning Hamiltonian used here, and those deduced from 
the further resummed, Kerr-like EOB Hamiltonian proposed in 5]. In view of the comparison showed in Fig. |21 we expect that the two 
models are very close to each other. 
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In the upper panels of Fig. ^] we plot the accumulative energy loss SEh starting from /gw = 40 Hz as a function of 
the instantaneous GW frequency /, for (15 + 15)Mq (upper left panel) and (15 + 5)Mq (upper right panel) binaries, 
each for 4 sets of initial spin orientations with maximal spins: aligned (dash-dot-dot curves), antialigned (dotted 
curves), (^5^ , ; , (/is^ ) — (60°, 90°; 60°, 0°) (denoted by generic-up, dash-dot curves), and (6'si , 0Si ; ; ) ^ 
(120°, 90°; 120°, 0°) (denoted by generic- down, dashed curves), as well as for the non-spinning configuration (solid 
curves). We use both SEHPF(3,3.5) (dark curves) and SEHPF(2,2.5) (light curves) models. In each panel, we also 
use vertical grid lines to mark LSSO frequencies. For the SEIIPF(2,2.5) model, all our evolutions go beyond their 
corresponding LSSO frequencies. The situation is a bit different for the SEIIPF(3,3.5). Indeed, as was shown in 
and in Fig. O above, the 3PN-E0B LSSO for mostly aligned fast-spinning BH's is drastically drawn inwards 
towards very high orbital frequencies. So high, indeed, that, for aligned and generic-up configurations, they fall out 
of the frequency range plotted in Fig. 1121 As a consequence, for aligned and generic-up configurations, the dynamical 
evolutions become rather non-adiabatic even before the formal LSSOs is reached. 

As we see from the plots, within a given GW frequency interval, binaries tend to emit more energy in configurations 
where spins are more aligned with the orbital angular momentum. This agrees with the results of and of Fig. 
above, showing that more aligned configurations are drawn towards more deeply bound states. This, together with 
the fact that LSSO frequencies are pushed higher in aligned configurations (as we also see from Fig. 0, can make 
the total energy releases in aligned configurations several times more than those in anti-aligned configurations. In 
Table IVIll we list values of SEh/M, accumulated from 40 Hz up to LSSO frequency (if reached) or ending frequency, 
otherwise, for configurations plotted in Fig. ^1 We also list the energy released below 40 Hz, and values of energy 
release when 2PN Hamiltonian and 2.5PN Pade flux are used. For (15 -I- 15)M0, maximally spinning binaries, the 
energy released from 40 Hz up to the end of our evolution (determined by one of the criteria (|4.1a|) - (|4.1d|) ') can range 
from 0.6% of M in the (anti-aligned configuration) to 5% of M (anti-aligned configuration), with the non-spinning 
configuration releasing 1.6% - 1.8% of M (in which 0.8% - 1.1% of M is released before LSSO). For (15 + 5)Mq 
binaries, the range is similar, from 0.5% to 5% of M, with non-spinning configuration releasing 1.2% ~ 1.4% of M (in 
which 0.7% - 0.8% of M from before the LSSO). We also note that the energies of around 0.8% of M and 0.5% of 
M are released below 40 Hz, for (15 -I- 15)M0 and (15 -I- 5)Mq binaries, respectively. See also Eq. (4.1) of ^ for an 
approximate analytical estimate of the energy released down to the LSSO, as a function of both rj — 17111112/ M'^ and 

XL- 

A. Evolution of the dimensionless rotation parameter J/E^ 

An important consistency check of the FOB approach concerns the dimensionless total angular momentum ratio 
|J|/i5^ ultimately reached by spinning black hole binaries. Indeed, if the FOB method would, at the end of its vahdity 
domain, predict a ratio |J|/i?^ larger than unity, this would preclude to match this end state with the newly born 
Kerr black hole expected from the coalescence of the two initial (spinning) black holes. This issue was investigated 
in the adiabatic approximation in There, it was shown that, when using the 3PN-accurate FOB Hamiltonian, the 
ratio \3\/E'^ estimated at the LSSO, was always smaller than unity. This result was not at all guaranteed in advance, 
and resulted from a delicate competition between the linear increase of |J| when increasing the (aligned) spin of the 
individual BH's, and its non-linear decrease because of the displacement of the LSSO towards smaller radii for aligned 
spinning configurations. It was found in |5i that the maximum of (| J|/i?^)LSSO was about 0.83, and was reached for 
XL ^ +0.3. 

The fact that this maximum value is significantly below one, leaves room for not running into any consistency 
problem even when taking into account the further changes of both E and J during the plunge that follows the 
crossing of the LSSO. Though our present attack on the problem does not properly consider the final matching 
between the plunge and the formation of a final Kerr hole, it goes beyond the previous treatments in going beyond 
the adiabatic approximation. In the lower panels of Fig. lT^ we plot the continuous time evolution of the ratio \3\/E'^ 
during the late stages of the inspiral and its subsequent non-adiabatic ending (which, in many cases, except in fact for 
the most dangerous aligned, is a post-LSSO plunge). It is convenient to use the gravitational wave frequency /gw to 
label the "time" along this evolution. Satisfactorily, we observe that, for all the binaries we have considered, |J|/i?^ 
decrease to below 1 before either the LSSO or the end of our evolution, whichever comes first. This means there 
are no a priori obstacles to having a Kerr black hole form right after the end of the non-adiabatic "quasi-plunge" . 
This means also that, contrary to an early suggestion [E^ based on rather coarse estimates, there is no ground for 
expecting a large emission of gravitational waves between the plunge and the merger. In Table IVIll we also list the 
values of \J\/E'^, at LSSO frequency (if reached) and ending frequency, for SFP(2,2.5) and SEP(3,3.5) models. 

In a pioneering work. Baker et al. [28j evaluated by a 3D numerical simulation the energy radiated from moderately 
spinning BH binaries with spins aligned or anti-aligned with the orbital angular-momentum. They started the (very 
short) numerical evolution close to the LSSO predicted by the effective potential method of Pfeiffer et al. As 
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already mentioned above, the numerical initial data chosen in these works rely only on an Initial Value Problem 
(IV P) formulation, and significantly differ both from the numerical initial data constructed by the HKV method 
|24Ll4a and from the predictions made by the EOB method (while the HKV and EOB results are quite close to each 
other I23L li^l. For instance. 23] estimated that the ratio between the orbital frequencies cjivp/'j-'eob was about 2. 
One should probably wait until HKV-type initial data for spinning BH's are evolved until coalescence to meaningfully 
compare their results with the results derived above for energy releases within the EOB approach. However, to have an 
idea of the current distance between analytical estimates and numerical ones, we have determined the energy released 
between the LSSO and the final frequency, at 2PN and 3PN order, for two of the spin configurations investigated by 
Baker et al. For spins aligned (anti-ahgned) and Xi — X2 — 0.17 (xi — X2 = 0.25), we find that the energy released is 
~ (0.6 - 0.9)%M [~ (1 - 3)%M]. Baker et al. found - (1.7 - 1.9)%A/ and ~ (1.9 - 2.1)%M, respectively. It should 
be noted that the energy released evaluated by Baker et al. includes also the ring-down phase. Ours does not. 

B. Complete waveforms describing the non-adiabatic inspiral and coalescence of precessing binary black 

holes 

Finally, though we have not yet carefully studied at which stage we could meaningfully join our "quasi-plunge" 
evolution to the formation of a ringing BH, we have decided, to show the promise of a purely analytical EOB-base 
approach to follow Ref. in matching (by requiring first-order continuity of the emitted waveform) the end of 
our waveform (here defined by the first violation of the "adiabatic criteria" (|4.1a|l -- (|4.1bl) ') to a ringdown waveform 
generated from the lowest I = m = 2 quasi-normal mode of a Kerr black hole. We determine the mass and spin 
parameters of the final hole by the energy and angular momentum of the binary at the end of our evolution: 

A^BH = £^fi„ , OBH - [|J|/^^']fi„ ■ (6.3) 

In Fig. 1131 we plot the complete waveforms, so obtained, for non-spinning, and maximally spinning (15 + 15)AfQ 
binaries in the generic-up and generic-down configurations (these refer to the configuration at /gw — 40 Hz, the 
starting point of evolution). We have shifted these waveforms in time so that the end of inspiral evolutions all happen 
at t = 0. Notice that at this stage the waveform which includes the ring-down phase should be considered as an 
example. Indeed, by restricting ourselves to the quasi-normal mode I ~ m = 2, we have tacitly assumed that the total 
angular momentum at the time the ring-down phase starts is dominated by the orbital angular momentum. However, 
this is not generally the case when spins are present and the quasi-normal modes with / 7^ 2 might be excited, as well, 
and contribute to the waveform. A more thourough analysis is left for the future. 

VII. CONCLUSIONS 

We provided a first attack on the problem of analytically determining the gravitational waveforms emitted during the 
last stages of dynamical evolution of precessing binaries of spinning black holes, i.e. during the non-adiabatic ending 
of the inspiral phase, and its transition to a plunge. We reviewed the various available Hamiltonian descriptions of the 
dynamics of spinning black hole (BH) binaries, and studied (following 5]) the characteristics of the stable spherical 
orbits that exist when spin-spin effects are neglected compared to spin-orbit ones. We derived the contribution to 
radiation reaction (for quasi-circular orbits) which is linear in the spins. Our results agree with the corresponding 
recent results of p34| . We then used this analytical description of the radiation-reaction-driven inspiral of spinning 
binaries to construct non-adiabatic models of coalescing binary waveforms. We compared the various models and 
concluded, in confirmation of previous results, that our current "best bet" for a non-adiabatic model describing 
the transition from adiabatic inspiral to plunge is obtained by combining: (i) an effective one body (EOB) [y, 
0, 3PN-accurate |2li resummed Hamiltonian, including spin-dependent interactions with (ii) Padc-resummed 
radiation reaction force (including spin-terms). Conclusion (i) is rather robust, since as Fig.n]shows, the PN-expanded 
Hamiltonian does not show any LSSO and differs significantly from the PN-expanded analytically computed function 
E{n); conclusion (ii) is based on the assumption that the flux function in the equal-mass case is a smooth deformation 
of the test-mass limit result. Since in the latter case Pade approximants were shown 19, 55j to have better agreement 
with exact numerical flux functions, we would conclude that this is also true in the equal-mass case. 

Our main results, obtained by means of this "best bet" EOB model are: 

(1) An estimate of the energy and angular momentum released by the binary system during its last stages of 
evolution: inspiral, transition from inspiral to plunge, and plunge; 

(2) The finding (which confirms the conclusions of 5j) that the dimensionless rotation parameter j / E'^ is always 
smaller than unity at the end of the inspiral; 
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(3) The construction of complete waveforms, approximately describing the entire gravitational- wave emission process 
from precessing binaries of spinning black holes: adiabatic inspiral, non-adiabatic transition between inspiral and 
plunge, plunge, merger and ringdown. Following these waveforms were constructed by matching a quasi-normal- 
mode ringdown to the end of the plunge signal. These tentative complete waveforms are preliminary because we did 
not include here a careful study of how to join, in a physically motivated manner, the last stages of the plunge to the 
merger phase. They extend, however, the (better justified) complete waveforms constructed in 0] to the more genral 
case of spinning and precessing binaries. 
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(^Sl, (l>S2, Osi, 032) 


Mn^ (Hz) 


/lsso(Hz) 


Rgr/M /In"" (Hz) 


/l"s^'o(Hz) 


R-fln ^ / M Pmax,2 


pmax,5 


(10 + 10) A/0 


no spin 


289 


285 


4.8 


287 


285 


4.2 


0.9150 




(0°,0°,0°,0°) 


745 


2145 


2.5 


466 


2145 


2.8 


0.3750 




(180°,0°,0°,0°) 


290 


285 


4.8 


285 


285 


4.2 


0.9166 




(180°, 0°, 180°, 0°) 


145 


145 


8.1 


145 


145 


7.4 


0.5587 




(60°, 90°, 60°, 0°) 


633 


873 


2.6 


502 


502 


2.6 


0.4851 


0.5883 


(120°, 90°, 60°, 0°) 


280 


280 


4.9 


269 


269 


4.4 


0.5420 


0.9472 


(120°, 90°, 120°, 0°) 


187 


187 


6.7 


186 


186 


6.0 


0.6096 


0.9536 


(15 + 15)Mq 


no spin 


192 


190 


4.8 


192 


190 


4.2 


0.8137 




(0°,0°,0°,0°) 


497 


1430 


2.6 


311 


1430 


2.6 


0.4550 




(180°,0°,0°,0°) 


193 


190 


4.8 


191 


190 


4.2 


0.8148 




(180°, 0°, 180°, 0°) 


98 


97 


8.0 


97 


97 


7.4 


0.6403 




(60°, 90°, 60°, 0°) 


429 


758 


2.6 


347 


347 


2.4 


0.5067 


0.6290 


(120°, 90°, 60°, 0°) 


185 


185 


4.9 


180 


180 


4.4 


0.5112 


0.9456 


(120°, 90°, 120°, 0°) 


124 


124 


6.7 


124 


124 


6.7 


0.6868 


0.9684 


(15 + 5)Mq 


no spin 


267 


265 


5.2 


265 


265 


4.2 


0.6023 




(0°,0°,0°,0°) 


862 


1442 


2.3 


479 


1442 


2.5 


0.3268 




(180°,0°,0°,0°) 


176 


175 


7.0 


177 


175 


6.2 


0.5188 




(180°, O'', 180°, 0°) 


141 


140 


8.3 


140 


140 


7.5 


0.4445 




(60°, 90°, 60°, 0°) 


715 


796 


2.4 


425 


743 


2.4 


0.4478 


0.6111 


(120°, 90°, 60°, 0°) 


208 


207 


6.2 


208 


208 


5.3 


0.5471 


0.7496 


(120°, 90°, 120°, 0°) 


224 


224 


5.9 


225 


225 


4.9 


0.5735 


0.8360 



TABLE I: We list the overlaps between STHTF(3,3.5), used as target model, and SEHPF(3,3.5), used as template model, for 
several binary masses and initial spin orientations. The two black holes are assumed to carry maximal spins Xi — X2 = 1, but 
for comparion we also show the results in absence of spins. The evolution starts at /in = 30 Hz. In the first three columns we 
list the ending frequency, the LSSO frequency and the BH radial separation at tan for the template model and in the second 
three columns we show the same quantities but for the target model. The last two columns contains the overlap maximized 
only over 2 extrinsic parameters Pmax,2 and maximed over 5 extrinsic parameters pmax,5, as described in the text. [Notice that 
for the spin configuration (60°, 90°, 60°, 0°) and masses (10 + 1O)M0 and (15 + 5)Mq , due to a pole in the Pade-approximant 
flux, we apply the Fade resummation only to the non-spinning terms in the flux.] 
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(^31,032,^31,032) 


/In " (Hz) / 


SEP /TT^N 

L330(nzj 


Rin7M /If ^ (Hz) 




Jn^/M pmax,2 Pmax,5 


(10 + 1O)M0 


no spin 


242 


237 


5.6 


237 


237 


4.5 


0.4691 


(0°,0°,0°,0°) 


628 


628 


2.9 


628 


628 


2.6 


0.3170 


(180°,0°,0°,0°) 


237 


237 


5.6 


237 


237 


4.5 


0.4681 


(180°, 0°, 180°, 0°) 


139 


139 


8.4 


139 


139 


7.6 


0.6433 


(60°, 90°, 60°, 0°) 


367 


367 


4.1 


342 


342 


3.5 


0.4197 0.4882 


(120°, 90°, 60°, 0°) 


234 


234 


5.7 


229 


229 


4.7 


0.4220 0.6015 


(120°, 90°, 120°, 0°) 


173 


173 


7.1 


172 


172 


6.3 


0.6681 0.9556 


(15 + 15)M0 


no spin 


158 


158 


5.6 


158 


158 


4.5 


0.4880 - 


(0°,0°,0°,0°) 


419 


419 


2.9 


419 


419 


2.6 


0.4044 


(180°,0°,0°,0°) 


158 


158 


5.6 


158 


158 


4.5 


0.4885 


(180°, 0°, 180°, 0°) 


93 


93 


8.4 


93 


93 


7.6 


0.7140 


(60°, 90°, 60°, 0°) 


240 


238 


4.2 


241 


241 


4.1 


0.4549 0.5186 


(120°, 90°, 60°, 0°) 


156 


155 


5.7 


152 


152 


4.7 


0.4827 0.6767 


(120°, 90°, 120°, 0°) 


115 


115 


7.1 


116 


116 


6.2 


0.7227 0.9442 


(15 + 5)Mo 


no spin 


232 


232 


5.7 


233 


232 


4.7 


0.6111 


(0°,0°,0°,0°) 


608 


608 


2.9 


608 


608 


2.6 


0.2695 


(180°,0°,0°,0°) 


167 


166 


7.3 


166 


166 


6.5 


0.8720 


(180°, 0", 180°, 0°) 


136 


136 


8.5 


136 


136 


7.7 


0.6743 


(60°, 90°, 60°, 0°) 


352 


352 


4.2 


367 


367 


3.2 


0.2696 0.4978 


(120°, 90°, 60°, 0°) 


192 


192 


6.6 


191 


191 


5.7 


0.6566 0.8173 


(120°, 90°, 120°, 0°) 


169 


169 


7.2 


167 


167 


6.5 


0.6207 0.8970 



TABLE II: Overlaps between STHTF(2,2), used as target model, and SEHPF(2,2.5), used as template model, for several binary 
masses and initial spin orientations. The two black holes are assumed to carry maximal spins Xi = X2 = 1, but for comparion 
we also show the results in absence of spins. The evolution starts at /n = 30 Hz. In the first three columns we list the ending 
frequency, the LSSO frequency and the BH radial separation at for the template model and in the second three columns we 
show the same quantities but for the target model. The last two columns contains the overlap maximized only over 2 extrinsic 
parameters pmax,2 and maximed over 5 extrinsic parameters pmax,5, as described in the text. 



(^31,032,031,032) 


/?SKF,nol'L 
/fin 


(Hz) /l336 (Hz) Rfl„ 


VM /ir(Hz) / 


SEP 

LSSO v^^zj 




(10 + 1O)M0 


(60°, 90°, 60°, 0°) 
(120°, 90°, 120°, 0°) 


633 
187 


872 
186 


2.7 
6.7 


660 
186 


1211 
186 


2.7 0.9860 
6.7 0.9953 


(15 + 5)M0 


(60°, 90°, 60°, 0°) 
(120°, 90°, 120°, 0°) 


743 
178 


767 
177 


2.3 
7.0 


564 
179 


564 
179 


2.9 0.9839 
6.9 0.9969 



TABLE III: Effect of radiation-reaction force along L over the binary evolution and waveforms by comparing SEHPF(3,3.5) 
with no Fl, used as target model, and SEHPF(3,3.5), used as template model, for several binary masses and initial spin 
orientations. The two black holes are assumed to carry maximal spins Xi = X2 = 1- The evolution starts at /n = 30 Hz. In 
the first three columns we list the ending frequency, the LSSO frequency and the BH radial separation at ten for the template 
model and in the second three columns we show the same quantities but for the target model. The last two columns contains 
the overlap maximized only over 2 extrinsic parameters pmax.2 and maximed over 5 extrinsic parameters pmax.s, as described 
in the text. [Notice that for the spin configuration (60°, 90°, 60°, 0°) and masses (10 + 1O)M0 and (15 + 5)M0, due to a pole 
in the Pade-approximant flux, we apply the Pade resummation only to the non-spinning terms in the flux.] 
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(^31,032,^31,032) 


/In"" (Hz) flMoiti^) Rtn'/M /ir(Hz) /Si'so(Hz) R|r/M p^.^,. 


(10 + 10)Mo 


(60°, 90°, 60°, 0°) 
(120°, 90°, 120°, 0°) 


632 872 2.6 616 1252 2.6 0.8566 
185 185 6.7 186 185 6.7 0.9762 


(15 + 5)M0 


(60°, 90°, 60°, 0°) 
(120°, 90°, 120°, 0°) 


743 767 2.3 661 772 2.5 0.8232 
178 177 7.0 179 178 6.9 0.9913 



TABLE IV: Effect of Pade and Taylor flux on the binary evolution and waveforms by comparing SEHTF(3,3.5), used as target 
model, and SEHPF(3,3.5), used as template model, for several binary masses and initial spin orientations. The two black holes 
are assumed to carry maximal spins Xi — X2 = 1- The evolution starts at /]„ = 30 Hz. In the first three columns we list the 
ending frequency, the LSSO frequency and the BH radial separation at tun for the template model and in the second three 
columns we show the same quantities but for the target model. The last two columns contains the overlap maximized only 
over 2 extrinsic parameters pmax,2 and max:imed over 5 extrinsic parameters pmax.s, as described in the text. [Notice that for 
the spin configuration (60°, 90°, 60°, 0°) and masses (10 + 1O)M0 and (15 + 5)Mq, due to a pole in the Pade-approximant flux, 
we apply the Pade resummation only to the non-spinning terms in the flux.] 



(^31,032,^31,032) 


/fin (^^J /lSSO (.J^^J ^fin /fin (^Z) jfL330 (^Z) Hfi„ /M Pniax,5 


(10 + 1O)M0 


(60°, 90°, 60°, 0°) 
(120°, 90°, 60°, 0°) 


650 1257 2.6 633 872 2.6 0.9959 
185 184 6.8 186 186 6.7 0.9988 


(15 + 5)M0 


(60°, 90°, 60°, 0°) 
(120°, 90°, 120°, 0°) 


702 766 2.3 743 767 2.3 0.9823 
178 178 7.0 178 177 7.0 0.9979 



TABLE V: Effect of quadrupole-monopole (QM) interaction on the binary evolution and waveforms by comparing SEHPF(3,3.5) 
with QM interaction, used as target model, and SEHPF(3,3.5) without QM terms, used as template model, for several binary 
masses and initial spin orientations. The two black holes arc assumed to carry maximal spins xi = X2 = 1- The evolution starts 
at fin = 30 Hz. In the first three columns we list the ending frequency, the LSSO frequency and the BH radial separation at ttm 
for the template model and in the second three columns we show the same quantities but for the target model. The last two 
columns contains the overlap maximized only over 2 extrinsic parameters pmax,2 and maximed over 5 extrinsic parameters pmax.s, 
as described in the text. [Notice that for the spin configuration (60°, 90°, 60°, 0°) and masses (10 + 1O)M0 and (15 + 5)M0, 
due to a pole in the Pade-approximant flux, we apply the Pade resummation only to the non-spinning terms in the flux.] 



(^31,032,^31,032) 


/•SEP,adiab 
/fin 


(Hz) /lsso (Hz) Kfi„ 


7m /ir(Hz) 


/ffs^o(Hz) 


Ri^/M pn,a.,5 


(10 + 10)Mo 


(60°, 90°, 60°, 0°) 


636 


1185 


2.6 


633 


872 


2.7 0.9666 


(120°, 90°, 120°, 0°) 


186 


185 


6.7 


186 


186 


6.7 0.9932 


(15 + 5)M0 


(60°, 90°, 60°, 0°) 


699 


827 


2.3 


743 


767 


2.3 0.9665 


(120°, 90°, 120°, 0°) 


177 


177 


7.0 


178 


177 


7.0 0.9914 



TABLE VI: Effect of assuming that spins evolve adiabatically. We compare SEHPF(3,3.5), used as target model, and 
SEHPF(3,3.5) obtained by averaging the spin couplings over an orbit, as template model, for several binary masses and 
initial spin orientations. The two black holes are assumed to carry maximal spins Xi = X2 = 1- The evolution starts at /in = 30 

Hz. In the first three columns we list the ending frequency, the LSSO frequency and the BH radial separation at for the 
template model and in the second three columns we show the same quantities but for the target model. The last two columns 
contains the overlap maximized only over 2 extrinsic parameters pmax,2 and maximed over 5 extrinsic parameters pmax,5, as 
described in the text. [Notice that for the spin configuration (60°, 90°, 60°, 0°) and masses (10 + 1O)M0 and (15 + 5)M0, due 
to a pole in the Pade-approximant flux, we apply the Pade resummation only to the non-spinning terms in the flux.] 
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(^Sl, 031,^2, 4>S2) 


[5EH]f<40Hz/M 


/lsso (Hz) [5EH]t%fB/M '\J\/E^'^^^^ 


/fin [SEHj^'^yM \J\/E''-^^ 


(15 + 15)M0, 3PN 


nospin 


0.0082 


190 0.0107 0.82 


325 0.0182 0.77 


(0°,0°,0°,0°) 
(180° .0°, 180° ,0°) 
(60°, 90°, 60° ,0°) 
(120°,90°,120°,0°) 


0.0086 
0.0077 
0.0084 
0.0079 


(1430) 

97 0.0033 0.51 
(760) 

123 0.0054 0.74 


474 0.0527 0.96 
194 0.0064 0.47 
440 0.0352 0.91 
242 0.0101 0.70 


(15 + 5)M0, 3PN 


nospin 


0.0048 


265 0.0084 0.62 


484 0.0141 0.58 


(0°,0°,0°,0°) 
(180° ,0°, 180° ,0°) 
(60° ,90° ,60° ,0°) 
(120°,90°,120°,0°) 


0.0049 
0.0046 
0.0049 
0.0047 


(1442) 

140 0.0034 0.14 

(793) 

177 0.0049 0.62 


819 0.0493 0.95 
289 0.0054 0.11 
719 0.0294 0.91 
351 0.0080 0.60 



TABLE VII: Energy released and the magnitude of angular momentum through the evolution (with spin-spin terms ignored). 
For non-spinning binaries, and four configurations of maximally spinning binaries, wc give the energy released below 40 Hz, 

from 40 Hz up to the LSSO, and from 40 Hz up to to the end of the evolution. [In some cases the evolution stops before LSSO 
can be reached.] We also show the corresponding values of \J\/E'^. Note that these results do not include the ring-down phase. 



26 



-0.02 
S~ -0-03 



2PN (no spins) 
2PN antialigned 
2PN aligned 
3PN (no spins) 
3PN antialigned 
3PN aligned 



0.02 0.04 0.06 0.08 



0.1 0.12 0.14 0.16 0.18 0.2 

MQ. 





-0.005 

-0.01 
-0.015 

-0.02 
-0.025 

-0.03 
-0.035 

-0.04 





1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 










' VV^' 






— 2PN (no spin) 






■ ■ ■ 2PN anti-aligned 






■ - 2PN aligned 






3PN (no spins) 






3PN anti-aligned 






2PN aligned 












0.02 0.04 0.06 O.C 



0.1 0.12 0.14 0.16 0.18 0.2 

MQ. 



FIG. 1: The energy for circular orbits as function of the frequency evaluated using the PN-expanded Hamiltonian (left panel) 
and the PN-expansion of the analytically computed function given by Eq. 12.331 (right panel) at various PN orders for maximal 
spins and equal mass binaries. 
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a -0.03 
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■ ■■ EOB 3PN antialigned 

■- EOB 2PN aligned \ 

"0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

MQ. 



FIG. 2: The energy for circular orbits as function of the frequency evaluated from the EOB Hamiltonian at various PN orders 
for maximal spins and equal mass binaries. 
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0.00 
-0.02 
-0.04 
m"-0.06 
-0.08 
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- 


\ * """c^ 


EOB 2PN (Damoui- [5]) 
EOB 3PN (Damoui- [5]) 

EOB 2PN (this paper) 
- ■ - ■ EOB 3PN (this paper) 

Adiabatic Taylor 2PN 
" - - Adiabatic Taylor 3PN 
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\i 
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5Cl 
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EOB 3PN (Damour [5]) 
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a 



^0.20: 

0.15^ 

0.10: 
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0.00 




-0.75 -0.5 -0.25 0.25 0.5 0.75 
5Cl 



FIG. 3: The energy (left panel) and the frequency (right panel) at the LSSO as function of xl/M'^ = ScR ■ L/K'P in the 
equal mass case for EOB Hamiltonian and PN-expanded analytically computed function E{yi) [see right panel of Fig. 0. The 
horizontal dashed line in the right panel marks the highest LSSO angular frequency for BBHs with total mass in the range 
1O-4OM0, assuming the LIGO freotiencv band 40 < /'nw < 240^2. 
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FIG. 4: Signal-to-noise ratio versus binary total mass at 100 Mpc for equal-mass binaries -with LSSO determined by the 
3PN. T?nT> TToTvi;i+^r,;o„ 




FIG. 5: Newton-normalized flux in the equal-mass case with both BH spins aligned (and maximal X — Xi — X2) with orbital 
angular momentum when T-approximants (left panel) and (upper-diagonal) P-approximants (right panel) are used. 




FIG. 6: Newton-normalized flux in the equal-mass case with both BH spins antiahgned (and maximal X = Xi = X2) with 
orbital angular momentum when T-approximants (left panel) and (upper-diagonal) P-approximants (right panel) are used. 
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FIG. 7: Comparison between the T- and (upper diagonal) P-approximant Newton-normalized flux in the equal mass case at 
3PN and 3.5PN order. 
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FIG. 8: Oscillations in r when spin-spin interactions are present, in (10 -I- 1Q)Mq (left panel) and (15 + 5)M0 (right panel) 
binaries. Dark curves show r /{ruo) as functions of /gw when both spin-orbit and spin-spin interactions are take into account, 
while light curves show the same quantity when only spin-orbit interactions are included. 
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FIG. 9: Relative orbital phase measured with respect to Stot x 
evolution at 40 Hz, with (Ssi , <;/>Si ; , ) = (60°, 90°; 60°, 0°), 



Lat. For maximally spinning (15 + 15)Mq binaries, we start 
and orbital phases i^orb = and 7r/2, and plot the difference 




FIG. 10: Comparison between waveforms from configurations with different initial orbital phases. For a (15 + 15)Af0 maximally 
spinning binary, we start evolution at 40 Hz, with {6si,4'Si',ds2j4'S2) = (60°, 90° ; 60°, 0°), and orbital phases (j>ovh ~ and n/2, 
and compare the waveforms detected with {F+, Fx ; 6, V?) ~ (1, 0; 7r/4, 0). 
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FIG. 11: Comparison between the three different prescriptions, SEj (dashed curves), SEh (dotted curves) and 5Eh (solid 
curves), for calculating energy losses. We use Newtonian-order radiation reaction in the left panel, and Fade at 3.5PN order in 
the right panel. We use the EOB Hamiltonian at 3PN order 
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FIG. 12: Accumulative energy release (upper panels) and instantaneous values of |J|/i5^ (lower panels) of (15 + 15)Mq 
(left panels) and (15 + 5)A4q (right panels) binaries. LSSO frequencies for the anti-aligned, generic-down, and non-spinning 
configurations are shown in vertical grid lines, while LSSOs of generic-up and aligned configurations are above the ranges of 
our plots. [Spin-spin terms are not included in these evolutions.] We use the SEHPF(3,3.5) model. 
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FIG. 13: Inspiral waveforms (which end at t = in our plot) matched to ring-down waveforms for non-spinning (light 
solid curve), and half-maximally spinning (15-1- 15) Mq binaries (left panel) and (15-1- 5)Mq binaries (right panel) in the 
generic-up (dark solid curve) and generic-down (dark dashed curve) configurations. We start our evolutions at 40 Hz, and use 
{F+,Fx;0,'p) = (1, 0; 7r/4, 0). In the plot we mark the position of the LSSO with solid curves. The waveforms have been 
shifted in time such that the end of the inspiral occurs at t/M = 0. 



